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^ ' 1. Introduction 



ABSTRACT 



Aims. We study timing properties of a large sample of y-ray bursts (GRB) detected by the Anti-Coincidence Shield (ACS) of the SPI 
spectrometer of INTEGRAL telescope. 

Methods. We identify GRB-like events in the SPI-ACS data. The data set under investigation is the history of count rate of the 
SPI-ACS detector recorded with a binning of 50 ms over the time span of ~ 10 yr. In spite of the fact that SPI-ACS does not have 
imaging capability, it provides high statistics signal for each GRB event, because of its large effective area. 

Results. We classify all isolated excesses in the SPI-ACS count rate into three types: short spikes produced by cosmic rays, GRBs and 
Solar flare induced events. We find some ~ 1500 GRB-like events in the 10 yr exposure. A significant fraction of the GRB-like events 
identified in SPI-ACS occur in coincidence with triggers of other gamma-ray telescopes and could be considered as confirmed GRBs. 
We study the distribution of durations of the GRBs detected by SPI-ACS and find that the peak of the distribution of long GRBs is 
at - 20 s, i.e. somewhat shorter than for the long GRBs detected by BATSE. Contrary to the BATSE observation, the population of 
short GRBs does not have any characteristic time scale. Instead, the distribution of durations extends as a powerlaw to the shortest 
time scale accessible for SPI-ACS, < 50 ms. We also find that a large fraction of long GRBs has a characteristic variability time scale 
of the order of 1 s. We discuss the possible origin of this time scale. 
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Gamma-ray bursts (GRBs) are the most violent explosions in the 
] universe. They are among th e first cosmic y-ray sour ces discov- 

■ ered more than 40 years ago (iKlebesadel et al.iri973h . In spite of 

■ the long observation history, the nature of the GRB phenomenon 
remains obscure. A number of missions were dedicated to the 

• study of the physical properties of GRBs. The isotropy of GRB 
distribution, indicating their extrag alactic origin, was est ablished 
by CGRO/BATSE observations dPaciesas et al.l Il999h. Long- 
lasting afterglows discovered by Beppo/SAX jPiro et al.lll998h 
confirmed the cosmological nature of GRB phenomenon and 
provided an insight in the energy budget of GRBs comparable to 
the energy release in the core collapse of massive stars. Complex 
time structure of GRB lightcurves, with variability dow n to mil- 
liseco nd time scale was revealed by CGRO/BATSE (Bh at et aP 
|1992|) . Existence of an additional delayed emission component 
in the GeV energy band was discovered by CGRO/EGRET 
( Schneid et al.lll992h and was recently confirmed by Fermi/LAT 
JAbdo et al.ll2009l) . 

The transient nature of GRBs makes the information avail- 
able about each individual event usually very fragmentary: a 
complete picture of the GRB phenomenon could not be estab- 
lished based on the observations of individual events. Instead, 
the complete picture could be obtained via population study in 
which different events appear as manifestations of a single GRB 
phenomenon. The population study should in the first place clar- 
ify if all the observed GRBs are produced by one and the same 
physical phenomenon or if the set of GRBs consists of sev- 
eral sub-populations which might possibly have different phys- 
ical origin. This problem is not simple: the quality of the data 



on more numerous weak bursts is not sufficient for quantitative 
analysis, while the bright bursts are not frequent enough to form 
a good statistical sample. 

Up to now, the largest sample of GRBs available for 
a statisti cal study ha s been collected b y the CGRO/BATSE 
(Paciesas 'eranil999l) and Konus-Wind JAptekaretal] Il995t) . 
One of the main findings of BATSE and Konus-Wind was the 
discovery of a dichotomy in the distribution of GRB durations 
(lMazetsetal.lll98U iKouveliotou et alJll993h . This distribution 
reveals two broad bumps. The first bump is at the durations 
Tgo < 2 s (where Tgo is the time over which 90% of the GRB 
fluence is registered), corresponding to the "short" GRBs. The 
second bump is at Tgo > 2 s, corresponding to the "long" GRBs. 
In the BATSE GRB sample, the population of long GRBs peaked 
at Tgo ^ 50 s, w hile the typ ical time scale of the short GRBs 
was Tgo ^ 0.3 s jKouveUotou et al. 1993) . Further analysis of 
a larger sample of BATSE bursts using improved duration esti- 
mates found that the durations of the long bursts s eem to peak in 
30 - 40 seconds range rather then at 50 seconds (iPaciesas et al.l 
[T999I) 

The division of GRB in two sub-populations is supported by 
observation of differences in spectral properti es of the two sub- 
classes: the short GRBs have harder spectra (IKouveliotou et al.l 
1993 [). Long GRB s are characterized by larger spectral lags 
dNorris et al.ll200ll) . The spatial locations of the bursts from the 
two sub-classes in their host galaxies are also different: long 
GRBs are tracing the star formation while the short GRBs are 
not correlated with the star formation regions (iBarthelmv et al.l 
20(31). 

The short versus long dichotomy is commonly believed 
to be related to the difference in burst progenitors. The long 
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GRBs seem to be associated with the final stage of evolu- 
tion of massive stars: gravitational collapse of the stellar core 
and/or supernova explosion (MacFadven & Woosley 1999j). The 
short GRBs are possibly related to coa lescence of compact ob - 
jects - neutron stars or black holes (iBarthelmv et al.l l2005b . 
In a number of cases sup emovae associated with long GRBs 
(IWooslev & Bloomll2006h were directly observed. Direct con- 
firmation of the hypothesis of association of the short GRBs to 
the merger events is still missing. 

Some problems of the division of the entire GRB populaiton 
in two classes remain ( Zhang 20111) . In particular, tight upper 
limits on the optical flux of several long GRBs exclude super- 
nova a s the GRB progen itor (iFvnbo et al.ll2006t iGal-Yam et al.l 
l2006t IValleetal.ll2006l) . A pecuhar long GR B with firmly 
exclu ded supernova progenitor, GRB 0606 14A dGehrels et al.' 
I2OO6I) has also negligible spectral lag typical for the short bursts. 
It is also possible that a separate population of softer spectra 
GRB-like events exists, the X-ray flashes (XRF). It is, however, 
not clear if the XR Fs are just a softer tail of the GRB distribution 
(Heise et al. 2001) or whether they are a different phenomenon 
tSoderberg et al...2006i) . Large fraction of short GRBs may con- 
tain softer extended emission comp onent, making their total du- 
ration as large as tens of seconds (iNorris & Bonnellll2006l) . In 
such situation the possibility of division of the GRB popula- 
tion in "sho rt" and "lon g " sub- classes becomes problematic. As 
claimed by ICline et aP (l2007h . the shortest bursts, with dura- 
tions of less than 100 ms, may also posses remarkably distinct 
properties: unusual spectra and non-isotropic angular distribu- 
tion (other studies of anisotropy of GRB subpopulations did not, 
however, come to the same conclusion (iBriggs et al.lll996l) ). a 
luminosity distribution suggesting Galactic origin. A popula- 
tion of Galactic sources, soft gamma- repeaters (SGR), is known 
to produce sh ort GRB-like events (iHurlev et al.l l2005l l2009t 
ISavchenko et al. 2010). The brightest of these events could be 
detected by existing instruments even if they occur outside the 
Milky Way, in the nearby galaxies. This type of events could 
constitu te a non-negligible fracti on of the observed short GRB 
sample (JHurlev et al.ll2005Ll2009l) . 

The focus of the present study is on the timing proper- 
ties of GRBs detected by the Anti-Coincidence Shield (ACS) 
of spectrometer SPI ori board of the INTEGRAL sattelHte 
(Ivon Kienlin et alJl2003ah . Although SPI-ACS does not have an 
imaging capability, it has a large effective area (up to about 
1 m^). SPI-ACS readout does not rely on any on-board trigger 
criteria and instead continuously measures the event rate all over 
the INTEGRAL mission lifetime with a time binning of 50 ms. 

The first atte mpt to systemati cally study GRBs with SPI- 
ACS was done bv lRau et all (l2005h using about 800 days of data. 
This study considered 179 GRBs also d etected by other s atel- 
Htes, but visible in the SPI-ACS. UnUke JRau et al.l (l2005h . we 
primarly rely on identification of the GRBs based solely on the 
ACS data, so that our GRB sample is not dependent on the data 
of other satellites. This allows us to perform a study of the timing 
properties of GRBs which is not affected by the selection biases 
(sensitivity, trigger configurations) of other telescopes. We sep- 
arately discuss a sub-sample of GRBs confirmed by the other 
satelUtes (Konus/Wind, Fermi/GBM, RHESSI, SuzakuAVAM, 
Swift/BAT, INTEGRAL/ISGRI, AGILE, IPN network). 

In the following sections we describe in detail the algorithm 
of detection of bursts in the SPI-ACS data and the classification 
of events into three classes: GRBs, instrumental short spikes and 
Solar flares. For each GRB we measure its duration and investi- 
gate if the GRB lightcurve possesses a characteristic time scale 



other than the overall duration. Based on this study we build the 
distribution of GRB durations and characteristic time scales. 

The main obstacle for the study of GRBs with SPI-ACS 
is the presence of the so-called "spikes" which have ligth- 
curves similar to the short GRBs, but are, mo st probably, asso - 
ciated to cosmic rays hitting the ACS detector (iRau et al.ll2005h . 
Association of the spikes with the cosmic ray events was ques- 
tioned by Del uit et al. ( 2004a) , who argued that the cosmic ray 
interactions in the SPI-ACS could not m anifest themselves as 
events with high count statistics. Instead, JDeluit et aP (l2004al) 
have advanced a hypothesis of possible existenc e of real very 
short G RBs missed by the other instruments. Later. ^Mina ev et all 
(I2OIOI) stacked the profiles of a sample of unidentified short 
bursts in the SPI-ACS and found a weak but recognizable af- 
terglow of the duration of more than a second. They interpreted 
this observation as an evidence in favor of the hypothesis of the 
association of the (fraction of) spikes with the real short GRBs. 

Our study of the origin of the sho rt spikes reported below, 
agrees with the principal conclusion of JRau et aP (l2005l) that the 
spikes are events induced by high-energy cosmic rays. At the 
same time, we show that the instrument effects of the ACS, the 
decay of cosmic ray induced radioactivity in BGO crystals, pro- 
duces the "afterglow" of the spikes discussed by iMinaev et al.l 
(l2010h . The detailed knowledge of properties of the spikes al- 
lows us to fully characterize them and to separate them from the 
real GRBs in our analysis. 



2. The INTEGRAL/SPI-ACS 



The spectrometer 
INTEGRAL 



SPI 



(IVedrenne et al.l l2003b onboard 
I Winkler et al.l (|2003|) is surrounded by an ac- 
tive ACS (Ivon Kienlin et al.ll2003bl) . consisting of 91 (only 89 
currently functional) BGO (Bismuth Germanate, BiiiGeT,0\2) 
scintillator crystals. The SPI-ACS also serves as a large effective 
area (up to ~ 1 m^) y-ray detector with a quasi-omnidirectional 
field of view dvon Kienlin et alj|2003al) . The design of the ACS 
is such that it has (almost) no sensitivity to the direction of 
y-ray s. 

The exact magnitude of the SPI-ACS effective area and its 
dependency on the direction and the energy is not well-known. 
It can be investigated through detailed simulations of the pho- 
ton propagation i n the d etector, as it was done, for example by 
iMereghetti et al.' ('2009"). However, this requires a mass model 
of the entire INTEGRAL satellite. An alternative method con- 
sists in making use of the events detected simultaneously by 
SPI-ACS and othe r detectors. This approach was exploited by 
IVigano & Mereghe tti (2009). 

The ACS data are event rates integrated over all the scintilla- 
tor crystals with a time resolution of 50 ms. The typical number 
of counts per 50 ms time ranges from about 3000 to 6000 (or 
even more during high Solar activity). 

The unique property of the SPI-ACS data is that, contrary to 
other existing GRB detectors, the readout does not rely on any 
trigger, so that a complete history of the detector count rate over 
the 10 yr period is recorded. This opens the possibility of an off- 
line search of GRBs in the compete data sample. 

The main challenge for the detection of GRBs in the coun- 
trate data of SPI-ACS is in the separation of real GRBs from 
other sources of variability of the signal, like the high-energy 
cosmic ray hits or charged particle events from the Solar flares. 
In the imaging detectors of GRBs, such separation could be done 
based on the directionality of y-rays contributing to the event. In 
the SPI-ACS separation of GRBs has to be done based on differ- 
ent selection criteria. 
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The burst like events are detected on top of a steady back- 
ground. The nearly constant background evolves only by a fac- 
tor of two over the ten years of operation. The background rate 
increases during the periods of passage of INTEGRAL sattellite 
through the Earth radiation belt and during the SPI "annealing" 
phases. These periods are excluded from the analysis reported 
below. 

The SPTACS data are searched for bursts in real time as 
part of the telescope operatio ns, using a dedicated IN TEGRAL 
Burst Alert System (IBAS) dMereghetti et al.ll2003h . The SPI- 
ACS triggers are u sed in the Inter-Planetary Network, the IPN 
dHurlev et al.ll2004b to derive the burst positions by the triangu- 
lation if other satellites of the IPN also detect the burst. 

Comparison of SPTACS to other similar scintillator type 
GRB detectors is given in Table [T] From this table one could 
see that the only GRB detector with the effective area compara- 
ble to that of the SPI-ACS was BATSE on board of the CGRO 



3. Detection of GRBs 

We use the SPI-ACS data of 1030 revolutions of the INTEGRAL 
satelUtes, from 20 to 1100 ((3230 days from 2002-12-12 to 
2011-10-16)). Some revolutions are excluded entirely because 
of the high solar activity or SPI annealing. The remaining total 
exposure time is 2412 days. 



3.1. Detection algorithm 

To identify burst-like events on top to the average background of 
SPTACS, we build the running mean of the ACS count rate on 
15 different time scales, from 100 ms to 10"* s. At each time scale 
6ts, the estimate of background count rate is obtained from the 
current value of the count rate on larger time scales, 6tb > St^. 
The burst-like events in the lightcurve with the time binning 5f , 
are identified as excesses in the count rate over the background 
count rate estimated in this way. Significance of detection of the 
burst detection depends, in general on the time scale 5tb used 
for the background rate estimate. We consider a burst-like event 
appearing on a time scale 6ts < 20 s as significantly detected if 
there is at least one scale Sth for which the significance of the 
signal detection on top of the estimated background exceeds a 
pre-defined threshold. The pair 5ts, 5ti, for which the significance 
of detection of the burst is the highest defines the "best time 
scales" and this highest detection significance is considered as 
the burst detection significance. 

This burst detection method includes an optimization for the 
signal/background time scales combination. This introduces a 
trial factor, which has to be taken into account in the significance 
estimate. Taking this into account, we set the detection signifi- 
cance threshold at a relatively high value, 6cr which corresponds 
to a; 5cr after the account of the trial factor 

The brightest bursts affect the running-mean estimates of the 
background even on the largest time scales. To correct for this 
effect, we perform the detection procedure twice, excluding the 
bursts detected at the first stage from the background estimation 
at the second stage. 

The background count rate is additionally affected by Solar 
activity. We evaluate the Solar activity related variability of the 
background via monitoring of the excess variance of the back- 
ground count rate on the "best" background time scale dth for 
every burst. The significance of the burst detection is corrected 
taking into account possible excess variance of the background. 



3.2. Duration measurements 

The most common and straightforward definition of the burst 
duration is the length of the time interval in which the burst sig- 
nal is detected as an excess over the background. With this def- 
inition, the measurement of the burst duration is obtained via a 
measurement of the significance of the burst signal in each bin 
of the lightcurve. The duration of the burst depends on the detec- 
tion significance threshold adopted in the analysis. This means 
that the burst duration defined in this way is not necessarily the 
characteristics of the burst itself, but depends on the sensitiv- 
ity of the instrument. Another drawback of this definition of the 
burst duration is that it could not take into account correlations 
between adjacent time bins. This leads to systematic underesti- 
mation of the duration: the burst could be not detectable in the 
narrow time bins, but still be significantly detectable on several 
time bins scale. 

An alternative techniq ue used to est imate burst durations is 
Bayesian Block analysis (IScarglelll998h . It consists in approx- 
imating the burst lightcurve with a sequence of intervals with 
constant flux. At first one interval is assumed. In each step of 
the iterative procedure one boundary is introduced. The itera- 
tions stop, when no more divisions are supported by the likeli- 
hood ratio test. This method has the advantage of not relying on 
an a-priori known background. This is particularly useful when 
the background is rapidly changing like in BAT and GBM. One 
of the disadvantages of the method is that it adopts a model of 
the signal in each block (constant) which is not physically mo- 
tivated. This could lead to an underestimation of the duration in 
the case of long and weak bursts. 

The method in which the burst duration is calculated based 
on significance of detection in fixed time bins could be improved 
to take into account weaker signals in several adjacent time bins. 
For this we run the analysis on several different pre-defined time 
scales. The set of the time scales is the same as used in the de- 
tection, see Section 13.11 We calculate the error on the burst du- 
rations from the comparison of the durations obtained assuming 
different detection thresholds, (4, 3 and 2 sigma). In general, du- 
rations calculated using different binning of the lightcurve are 
different. We compare the durations calculated for different bin- 
nings and choose the time scale which provides the minimal un- 
certainty to measure the start and stop of the burst. 

Similarly to the bayesian blocks methodology, the methods 
adopting binning of the lightcurves implicitly assume constant 
rate within the bin. Physically motivated models of the bursts 
profiles can be more successful in interpreting the properties of 
the burst prompt emission light curve. For exampl e, profiles of 
the b urst may be fitted with a collection of pulses ( iNorris et al.l 
il996|). However, more complex models the compromise gener- 
ality of the analysis. 

To verify the precision of the measurement of the burst du- 
rations, we perform an analysis based on the fixed time bins and 
on the Bayesian blocks and compare the results. We find that 
the burst durations found using the two methods are consistent 
within the uncertainties. 

The duration methods discussed above are based on statis- 
tical testing of the signal over the background. Since different 
instruments have different sensitivity to the signal and different 
background levels, the durations measured in this way depend 
on the instrument performance. A definition of the burst dura- 
tion which is less dependent on the instrument sensitivity is the 
so-called Tyo-duration, which is equal to the time during which 
90% of the burst flux arrives. To measure Tgo we identified the 
total flux within the detected duration, and adjust the start and 
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Suzaku/HXD-II WAM 


CGRO/BATSE LAD 


Beppo-SAX/PDS 


GRBM 


INTEGRAL/SPI-ACS 


Fermi GBM 


Crystal 


BGO 


Nal(Tl) 


CsI(Na) 




BGO 


BGO 


Energy range (keV) 


50-5000 


20-2000 


40-700 




>75 


150-30000 


Effective area (cm2) 


800 @ 100 keV 
400 @ 1 MeV 


2000® 100 keV 
150@lMeV 


700 @ 200 keV 
100@ 1 MeV 




3000@100keV 
< 8000® IMeV 


120@200 keV 
120@lMeV 


Time resolution 


31.25 ms 


2 ms 


7.8 ms 




50 ms 


5A/S 



Table 1: Comparison of SPI-ACS to other GRB detectors. 



the stop until 5% of the flux is removed from the beginning and 
from the end of the burst. 



3.3. Types of burst-like events in SPI-ACS 

The total number of burst-like events found in our analysis is 
8223 1 . Most o f these events ar e a special type of events called 
"short spikes" dRau et al.ll2005h . At long time scales, variations 
of the SPI-ACS count rate are mostly due to the changes of Solar 
activity. Below we explain the algorithm used to separate the 
GRBs from the spikes at short time scales and from the Solar 
flares at long time scales. 



3.3.1. Solar activity 

The Solar activity in the SPI-ACS introduces variability of vari- 
ous magnitudes, ranging from the minor flickering and slow evo- 
lution to giant flares increasing the SPI-ACS background by fac- 
tor of few for days. As mentioned above, the variability of the 
background due to the solar activity poses a problem for the cor- 
rect calculation of significance of the detection of GRBs. Our 
procedure for the calculation of significance of the burst takes 
into account the excess variance of the background during the 
periods of Solar activity. 

Lightcurves of some solar flares could sometimes appear 
similar to those of GRBs. It is not possible to distinguish such 
flares from GRBs based on the SPI-ACS data alone. To distin- 
guish the Solar flares from the GRBs we use the data of moni- 
toring of the Solar activity by the solar weather prediciton center 
(SWPC) Q to exclude time intervals which might be affected by 
the Solar activity. 




Seconds since the peak 



Fig. 1: The stacked profile of the spikes (blue data points). The 
extended tail is fitted with the decay of residual cosmic ray in- 
duced radioactivity of the BGO found from FLUKA simulations 
(green curve). 

Alternatively, the spikes (or a fraction of them) could be real 
GRBs which escape detection by other instruments. This could 
arise because of large effective area of ACS at the energies above 
~ 1 MeV. If the spikes are real GRBs, they should have very hard 
intrinsic spectrum, to avoid detection by other instruments. The 
hard and very short bursts would be a phenomenon of particu- 
lar interest. A distinct population of bursts with similar proper- 
ties might have been found in the Swift/BAT, BATSE and Konus 
Wind burst samples dCline etal.i ri999. 2007). 

We dedicate the next Section to study of the Short Spike 
properties. 



3.4. Short spikes 

The majority of burst-like events in SPI-ACS are "short spikes", 
with typical durations of about the about 50 - 100 ms, i.e. one 
or two bins of the lightcu rve. The rate of these events is very 
high - about 30 spikes/day dRau et al.ll2dollDeluit erani2004bh . 
Spikes are usually detected only in one 50 ms bin and very in- 
frequently in two, which suggests that the real duration is much 
smaller than the bin size. 

The nature o f these events was controversial. The analysis by 
iRau et all d2005h suggests that the spikes are produced following 
hits of the ACS by high-energ y cosmic rays. The observed flux 
of high-energy cosmic rays is dAdriani et al.ll20I lb 



CR 



10^ 



-'CR 



lOTeV 



-1.7 



m ^sr 'd ' 



(1) 



Taking into account the effective area of ACS, the observed rate 
of spikes suggests that the events are produced by cosmic rays 
with energies above ~ 30 TeV. 



http ://w w w. swpc.noaa.gov/ 



4. The nature of the short spikes 

The properties of the short spikes reveal a remarkable stability 
over the whole mission span (10 years of data). While the back- 
ground is slowly changing with the solar cycle, the sensitivity- 
corrected spike rate remains constant. This observation makes 
a Solar origin of spikes unlikely, unless there is a high-energy 
solar phenomenum not modulated with solar cycles. If cosmic 
rays are resposible for the spikes, however, they should not be 
affected by the solar activity and thus have energy high er than 
> 10 GeV, where the Sol ar modulation effect vanishes dParkeJ 
[T95l[Shikaze et al.ll200l 

To characterize the time profile of the spikes we stack the 
lightcurves of all the very short bursts detected by SPI-ACS 
(shorted then 0.5 s) shifted in time to match the bin with the 
highest countrate. 

The stacking analysis reveals a clearly recognizable weak 
(0.1 % of the total number of counts or just a few counts in 
a single spike on average) afterglow lasting for almost 3 sec- 
onds (Figure[T]i. From this figure one could see that the "prompt" 
phase of the spike is characterized by an exponential rise / ex- 
ponential decay profile. The rise and decay time are constrained 
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to be shorter than ~ 25 ms, half of the time bin of the ACS 
Ughtcurve. The overall event number and count statistics in the 
stacked sample is largely dominated by spikes, so that the pres- 
ence of confirmed short GRBs or SGR flares in the sample does 
not affect the average spike profile found from the stacking anal- 
ysis. To verify that the long tail of the stacked profile is not due to 
the presence of a sub-population of rare real GRBs in the stacked 
sample, we studied distribution of countrates in individual events 
in the time bins of the tail. This distribution follows a Gaussian 
with the center displaced from zero. If the tail of the profile 
would be due to the presence of a small number of real GRBs, 
the distribution would be instead a Gaussian centered at zero, 
with a small excess at positive count rates. We have also verified 
that the profiles of individual brightest spikes follow closely the 
averaged stacked profile, with a clearly visible tail. 

iRau et all (l2005h suggested that the spikes are produced fol- 
lowing energy deposition in the ACS by a high-energy cosmic 
ray. The cosmic ray initiates a shower of secondary particles 
through pair production and bremmstralung and then propagat- 
ing further to the Ge detectors of SPI and saturating them. One 
of the arguments presented by iRau et all (12005 ) in favor of in- 
terpretating of the spikes as events produced by the high-energy 
cosmic rays is the correl ation of the spike s with events of satura- 
tion of the SPI detector iRau et all (l2005h noticed that saturation 
of at least one Ge detectors of the SPI often (but not always) 
accompanies a short spike and that saturation of more than one 
SPI detector is always accompanied by a spike-like increase in 
SPI-ACS count-rate. 

Our analysis of the correlation between the spikes and the 
sa turation of SPI de tectors does not fully confirm the observation 
of iRau et alj ( |2005|) . Although we find a significant correlation 
between the spikes and SPI saturation moments (1600 coinci- 
dence events with expected 40), we find that saturations of more 
than one SPI detector are not always accompanied by spikes. 

In particular, the saturation of all the SPI detectors is always 
accompanied by a gap in the SPI-ACS data as well, so that no 
spikes could be recorded. Although there is a significant cor- 
relation between spikes and "cosmic-ray track-like" saturation 
events in SPI (more than 3 detectors in a row are saturated), some 
60 % of the track-like events are not accompanied by spikes. The 
correlation is still weaker in the saturation events which do not 
have track-like appearance. Only 5% of those events correlate 
with the spikes. 

Problems of this hypothes is were highlighted bv lDeluit et al.l 
(l2004bl) . iDeluit et ail (i2004ah have noticed that all the particles 
from a cosmic ray induced shower in ACS would be counted 
as just one particle, because of the finite time resolution of the 
ACS counting electronics (0.6 /ys). Indeed, the shower particles 
travel with a speed close to the speed of light. The transit time 
of the ACS is much shorter than the time resolution of the de- 
tector This means that the cosmic ray induced shower could not 
produce a spike in the ACS count rate. The shower particles also 
could not be responsible for the second-long afterglow shown in 
Fig.[Tl as highhghted bv lMinaev et alj ( l2OI0l) . 

A potential way out of this problem is that the cosmic ray 
induced shower in ACS is able to excite long-lived phosphores- 
cence of the detector In this case the phosphorescence signal 
persists on longer time scale and could be time-resolved by the 
ACS front-end electronics. This effect is actually observed in an- 
other detector on board of INTEGRAL, PICsIT, in which strong 
(~ 2500 counts s"') and short (< 1 70 ms) spikes are as sociated 
with cosmic ray induced showers (ISegreto et al.l 120031) . Time- 
tagging of photons and energy resolution of PICSIT allow iden- 



tification of the spikes as the phosphorescence of the detectors 
following the cosmic ray hits. 

SPI-ACS uses undoped BGO scintillators, featuring very 
rapid deexcitation - with a characteristic times of about 60 - 
300 ns followed by weak, 0.005%, residual phosphorescence on 
the time scale of 3 ms. If the energy deposit by a cosmic ray is 
large enough, even weak residual phosphorescence signal could 
be above the trigger threshold of the ACS front-end electronics. 
In this case, the detector will continue to count events on the time 
scale up to the ~ ms following the CR hit. This would explain 
the appearance o f a spike in the ACS count rate as suggested by 
IRau et all ( l2005h . 

However, prolonged phosphorescence of BGO crystal could 
not be responsible for the spike afterglows shown in Fig. [T] be- 
cause persistence of the phosphorescence signal on the ~second 
time scale would produce a flat constantly high countrate R ~ 
50 ms/0.6 yus ^ 80000 counts/bin on this time scale, rather 
than "spike-with-a-tail" lightcurve. Instead, the tail could be 
attributed to another feature of the BGO crystal, the induced 
radioactivity. Activation of the BGO crystals with the cosmic 
ray protons is known to contribute to the measured background 
dTruscott & Dveijri997l) . We have investigated the radioactivity 
induced in the BGO crystals by high-er iergy cosmic ray pas- 
sage using FLUK A Monte-Carlo code jBattistoni et al.li2007l: 
iFerrari et al.l 120051) . These simulations reveal that the main ra- 
dioactive decay channel is the decay of ^Li with the half-life time 
838 ms. There is also a number of other isotopes with compara- 
ble decay times although lower abundances. A fit of the observed 
time profiles of the extended tails of the spikes with the time pro- 
files of decay of the cosmic ray induced radioactivity found from 
FLUKA simulations is shown in Fig.[T] One could see that the 
time profile of the tail is well explained by the exponential decay 
of the induced radioactivity in the BGO crystal. 

To summarize, the observed properties of the spikes are con- 
sistent with the assumption that they are produced by high- 
energy cosmic rays passing through the ACS and inducing phos- 
phorescence and radioactivity in the BGO crystals. 

Knowledge of the physical origin of the spikes enables us 
to efficiently reject them in the analysis of the GRB sample. All 
the spikes are expected to have identical time profiles, with the 
overall normalization determined by the energy of the primary 
cosmic ray. Rejecting all the events which have the time profiles 
consistent with that of the template spike profile shown in Fig.[T| 
will remove most of the spikes, while suppressing only a small 
fraction of the real GRBs which occasionally have time profiles 
similar to the spikes (a pronounced peak followed by a tail with 
fixed fractional amplitude and decay time). Taking this into ac- 
count, we compute the likelihood for each event to be a spike 
rather than real GRB. We calibrate this likelihood with simu- 
lations and choose the threshold to reject the spikes in the real 
data. 



5. Properties of GRBs detected by SPI-ACS 

5.1. GRB durations 

Rejection of spikes and Solar flare events from the overall sam- 
ple of burst-like events in SPI-ACS leaves a sample of 1416 GRB 
lightcurves available for the timing analysis. Figure |2]shows the 
distribution of durations Tgo for these GRBs. This distribution is 
compared to the distribution of durations of the spikes shown by 
the blue histogram in Fig.|3] 

One could see that the spikes dominate the total number of 
events with durations up to several seconds. However, the char- 
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Fig. 2: Distributions of the Tgo durations. Dark and light shad- 
ings of the filled region correspond to 5 and 3 sigma confidence 
level for the event not being a spike. Solid black curve is for the 
confirmed bursts. The distribution of durations of GRBs from 
BATSE 4B GRB sample, scaled to match the normalization of 
the ACS sample, is shown by the black dashed line. 



acteristic time profile of the spikes allows their efficient rejec- 
tion even on the shortest time scales. In particular, a significant 
number of GRBs with durations down to ~ 100 - 150 ms is de- 
tected in the SPI-ACS data. Such bursts are identified as events 
for which the lightcurve is not consistent with the spike time 
profile. Two shaded histograms in Fig.|2]show two different con- 
fidence levels for an events not to be spikes. One could see that 
reducing the acceptance threshold leads to a sharp increase of 
the number of GRB-like events on the shortest time scale where 
difference between the spikes and real GRBs is more difficult to 
characterize. Taking this into account, we restrict our attention 
to GRBs with durations longer than 100 ms in the following. 

Distribution of durations of the spikes follows a powerlaw on 
large time scales. This is explained by the fact that the span of the 
time interval during which the cosmic ray induced radioactivity 
produces a detectable signal in the BGO crystal scales propor- 
tionally to the overall number of the radioactive atoms which 
is in turn proportional to the energy of the primary cosmic ray. 
Taking into account that the spectrum of cosmic rays follows 
a powerlaw, one could find that the distribution of durations of 
the spikes should also follow a powerlaw. Amplitude of statis- 
tical fluctuations of the spike duration distribution determines 
an upper bound on the possible number of real GRBs in each 
duration time bin. This upper limit is shown in Figure |3] as an 
upper boundary of the pink-shaded region. Lower boundary of 
the pink shaded region is the distribution of the GRB-like events 
calculated with tight rejection threshold for the spikes, shown in 
Figure 12] 

Real GRBs are expected to produce signals not only in ACS, 
but also simultaneously in other GRB detectors. To verify the 
quality of our procedure of selection of GRB events from the 
set of burst-like events in SPI-ACS, we have searched for the 
events occurring in coincidence with triggers of other GRB de- 
tectors. Table|2]summarizes the result of our search. The number 
of events found in coincidence with triggers of other detectors is 
comparable to the number of events selected based on the SPI- 
ACS data alone. Slightly larger number of events found using 
SPI-ACS alone is readily explained by the large effetive area of 
ACS. In our search of events occurring in coincidence with trig- 
gers of other detectors, we have adjusted the search time window 



Fig. 3: Distribution of GRB durations (shaded region) compared 
to the distribution of durations of spikes (blue histogram). 



(within which the events are declared to be coincident) in such a 
way that the we expect 1 random coincidence. 

The black solid histogram in Figure |2] shows the distribution 
of durations of events occurring in coincidence with the triggers 
from other detectors. One could see that this distribution is con- 
sistent with that found based on the analysis of the SPI-ACS data 
alone, although with somewhat lower statistics. 

Although the presence of the short spikes reduces the sen- 
sitivity of the SPI-ACS to the shortest (< 100 ms) bursts, the 
number of short bursts with longer durations (100 ms - 2s) is 
rather well constrained. The fact that distribution of durations of 
the independently detected bursts is so similar to the distribution 
of the confirmed burst durations (Figure|2), shows that our spike 
filtering technique efficiently rejects most of the spikes and re- 
tains most of the bursts with durations longer than two ACS time 
bins. 



5.2. GRB peak fluxes 

To verify the efficiency of rejection of spikes based on the tem- 
plate profile, we study the distribution of peak countrates in the 
remaining GRB sample. The peak countrate measured by SPI- 
ACS provides an estimate of the peak flux of the GRB. SPI- 
ACS detector could not measure the peak flux of each detected 
GRB, because its effective area depends on the (in general, un- 
known) arrival direction of GRB photons. However, GRBs ar- 
rive isotropically, so that the average peak flux distribution could 
be measured for all the GRB population, using a known angle- 
averaged effecitve area of ACS. This distirbutiton is shown in 
Figure |4] One could see that the distribution has the form of a 
broken powerlaw with slopes -1.46 + 0.04 at high countrates 
and -1.12 + 0.01 at the countrates below 4000. This count 
rate corresponds to a GRB peak flux ~ 10 ph/cm-s in the en- 
ergy band above 80 keV or, equivalently, to an energy flux 
5 X 10"^ erg/cm^ s. For comparison, th e distribution of the GR B 
peak fluxes as measured by BATSE (iKommers et al.l [l998h is 
shown at the same figure. It is scaled for the convenience of 
the presentation. Notably, both the break and the slopes in the 
BATSE peak flux distribution are close to the one we observe. 



5.3. Characteristic time scale of long GRBs 

We have investigated the variability patterns of the GRB 
lightcurves, using conventional time series analysis tools, such 
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GRB 461 



452 292 145 127 21 

Table 2: SPI-ACS triggers detected simultaneously by other means 
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as power spectral density (PSD) and structure funciton analysis. 
These tools are well suited for the analysis of long GRBs, in 
which the dynamic range of time-scales / frequencies for each 
GRB event is large enough. 

We use the first order structure function to quantify the burst 
temporal properties. The structure function analysis is analogous 
to the autocorrelation function (ACF) analysis, successfully used 
by a number of authors in investigations of the G RB variability 
(lBorgonovoll2004t lin 'tZand & Fenimorel '1996'). We prefer the 
structure function analysis to the PSD analysis which is better 
suited for the study of variability of persistent, rather than tran- 
sient sources, because of the nature of the Fourier transform. 

The structure function analysis could be used to identify a 
characteristic variability time scale, which appears as a local- 
ized "deep" feature superimposed on the powerlaw type behav- 
ior. We use this property of the structure function to investigate 
the possibile existence of a characteristic variability time scale 
of GRBs. The existen ce of such a timesc a le in t he T„,r -Is 
range was claimed by iBeloborodov et all (l2000l) based on the 
PSD analysis of a set of BATSE GRBs. 

Our structurre funct ion analysis confirms the result of 
IBeloborodov et al](l2000l) . Figure|5]shows the characteristic time 
scales of long GRBs found form the structure function analy- 
sis and plotted as a function of the GRB duration. One could 
see that large fraction of GRBs has the characteristic time scale 

r„„, = 1 ± 1 s. 



6. Discussion 

Our analysis of the sample of GRBs detected by SPI-ACS shows 
that the distribution of GRB durations shown in Figure |2] has 
two components. Similarly to BATSE duration distribution, it 
has broad bump corresponding to the long GRB population. The 
typical duration of long GRBs in SPI-ACS (- 20 s) is somewhat 
shorter than in BATSE (^ 50 s). This is in agreeme nt with the 
previous studies of the SPI-ACS bursts performed by iRau et all 



Fig. 5: One the top panel: the first characteristic time scale vs du- 
ration (Tgo) for long GRBs. On the bottom panel; the distribution 
of the first characteristic time scales. 



( l2005h . The difference in the typical duration could be attributed 
to the difference in the energy ranges in which BATSE and 
SPI-ACS are sensitive. The dependence of the duration distribu- 
tion on the energy range was previously observed in Swift/BAT , 
BATSE, Beppo-SAX and HETE data dFenimore et al.l 1 1 9951: 
I Sakamoto et al .11201 ih . Longer bursts are known to have softer 
spectra. Besides, the spectra of individual GRBs have hard-to- 
soft evolution so that the burst duration grows with decreasing 
photon energy. 

Contrary to BATSE burst duration distribution, the distribu- 
tion of durations of SPI-ACS bursts, shown in Figures |2] and |3] 
does not reveal short bursts concentrated around a characteris- 
tic time scale. Instead, the burst duration distribution follows a 
powerlaw at short time scales Tgo < 3 s. The powerlaw extends 
to the shortest time scale accessible to the ACS (^ 100 ms) and, 
most probably, to the longer time scales, beyond ~ 3 s. 

This might be either due to the instrumental effects or might 
reflect the GRB physics. An instrumental effect could arise e.g. 
because of the difference of sensitivity between BATSE and 
SPI/ACS. If short GRBs typically have extended tails, which 
were not detected by BATSE, but are visible to SPI/ACS, part of 
the bursts which might have been responsible to the short burst 
duration distribution bump at 0.3 s, might appear in SPI/ACS as 
longer bursts in the ~ 3 s duration range. Otherwise, the differ- 
ence between the short and the long bursts may not be fully char- 
acterized by the criterion Tgo S 2 s, as suggested by BATSE dis- 
tribution. Separation in the hardness and spectral lags have to be 
used to establish the separation between the two main varieties 
of the GRBs. As discussed by Zhang ( 2011 ), it may be natural to 
adopt a different classification scheme: to Type I ("compact star 
GRBs") and Type II ("massive star GRBs") bursts. 



V.Savchenko, A.Neronov. TJ.-L.Courvoisier: Timing properties of gamma-ray bursts detected by SPI-ACS detector 



Although our analysis does not reveal a separate population 
of very short b ursts w ith the durations Tgo < 100 ms, discussed 
bylCline et al.l ( |2007|) . we could not completely rule out the exis- 
tence of such sub-population. Such very short bursts could be 
easily confused with the cosmic ray induced spikes in ACS. 
To distinguish real bursts from the cosmic -ray induced spikes 
a more advanced analysis based on the simultanenous detections 
by the other intruments, e.g. Suzaku/WAM, has to be done. 

The sample of the SPI-ACS bursts, presented in this work, is, 
of course, subject to observational biases related to the detection 
efficiency of our method. It was shown that in the case of BATSE 
specific detection efficiency issues have significant effect on the 
burst duration distribution, increasing the fraction of short bursts 
to about 40% dLee & Petrosianiri996h . Unfortunately, low num- 
ber of measurable burst characteristics accessible in SPI-ACS 
(limited to timing properties of lightcurves in a single energy 
band) make similar analysis impossible in the case of SPI-ACS. 
A detailed study of the SPI-ACS GRBs detected simultaneously 
with other satellites may help in recovering the intrinsic distri- 
butions of the GRB properties. Alternatively, selection biases 
might be assessed via a consideration of different trigger criteria. 
It has been shown that a sample of the short bursts may be bet- 
ter probed with the ti me-to-spill flux-bas ed rather the traditional 
fluence-based trigger (iNorris et al.llI984l ). 

The second result of our work is the identification of a char- 
acteristic time scale of variability of long GRBs, see Figure |5] 
This time scale can be roughly attributed to the typical separa- 
tion between the individual pulses or activity periods of GRBs. 
It is interesting to note that the characteristic time scale of long 
GRBs coinsides with the time scale at which a break in the distri- 
bution of GRB durations is observed (see Figure |3). This might 
indicate that many long GRBs are sequences of ~ 1 s pulses, 
with the overall burst duration determined by the overall number 
of pulses. The "pileup" in the distribution of GRB durations, at 
the time scale ~ 20 s shows that long GRBs typically consist of 
~ 10 - 20 pulses. 

The m ajority of GRB with me asured redshifts are situated at 
Z ~ 1 - 2 dSalvaterra et al.JEOl lb . This implies that the charac- 
teristic time scale in the GRB rest frame is of the order of 0.5 s. 

Our measurement of the characteristic pulse scale is ro- 
bust but simplistic. Extensive studies of the indi vidual pulse 
proper ti es were performed by a number of groups p^orris et all 
dT996l) : iHakkila & Nem iroff' ("200?); 'Hakkila &Preecd (l20IIh . 
In their pioneering work N orris et al. (1996), matched the ob- 
served lightcurves by a combination of model pulses. Notably, 
their result - the characteristic time scale of the order of a second, 
is close to the one we obtained with model-independed structure 
function analysis^ 

Results of ([Norris etaL 1 9961) and fu rther studies (e.g. 
IHakkila et al.ll2008l: iRvde & Svenssonll2002h also indicated that 
the pulse properties (including characteristic timescales) are cor- 
related. Since we use a wide range of the peak fluxes, it may 
should have an effect on our measurement of the characteristic 
1 second scale. However, this time scale is readily detectable in 
the SPI-ACS data even after smoothing introduced by ignoring 
the correlations. The robustnes of the detection supports the idea 
that it should be associated with a real physical characteristic of 
the bursts. 

The distinct characteristic time scale of the order of a second 
might be associated with the variability of the GRB central en- 
gine. In this case the remarkable persistency of the this time scale 
seem s to be naturally attributed to the similarity of the progen- 
itors (iBeloborodov et alll2000l) . Alternatively, it can be related 
to the size of the emission region. Further investigation of the 



nature of the characteristic time scale, exploiting the rich timing 
data of the SPI-ACS, is required to clarify its origin. 
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